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(54) Method for thermal printing 



(57) A method for printing an image using a thermal 
printing system comprising a thermal printer having a 
thermal head TH incorporating a plurality of energisable 
heating elements H n and a thermographic material m, 
comprises several steps. These steps comprise supply- 
ing image data Id corresponding to a plurality of output 
values D n to a processing unit, reading a first table LUT1 
comprising first entries for desired output values D k and 
second entries for sequences S k of time slices s ( of ac- 
tivation pulses, transformation of the image data Id into 
corresponding sequences S k of time slices of activation 



pulses, providing the time slices of activation pulses to 
the heating elements H n of the thermal head TH, and 
printing the image by transporting the thermographic 
material past and adjacent to the thermal head and by 
activating the heating elements H n of the thermal head. 
Herein, the transformation of the image data Id compris- 
es a step of establishing a relation between output val- 
ues D kcm or D kcm , which are both continuous and mon- 
otone and in the latter case linear as well, and corre- 
sponding sequences S k of time slices Sj of activation 
pulses. 
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Description 

FIELD OF THE INVENTION 

5 [0001] The present invention relates to a device operable for applying thermal energy to a recording medium, the 
device comprising a thermal head having energisable heating elements which are individually addressable. In partic- 
ular, the recording medium is a thermographic material, and the head relates to thermal imaging, generally called 
thermography. 

10 BACKGROUND OF THE INVENTION 

[0002] Thermal imaging or thermography is a recording process wherein images are generated by the use of im- 
agewise modulated thermal energy. Thermography is concerned with materials which are not photosensitive, but are 
sensitive to heat or thermosensitive and wherein imagewise applied heat is sufficient to bring about a visible change 

15 in a thermosensitive imaging material, by a chemical or a physical process which changes the optical density. 

[0003] Most of the direct thermographic recording materials are of the chemicat type. On heating to a certain con- 
version temperature, an irreversible chemical reaction takes place and a coloured image is produced. 
[0004] In direct thermal printing, the heating of the thermographic recording material may be originating from image 
signals which are converted to electric pulses and then through a driver circuit selectively transferred to a thermal print 

20 head. The thermal print head consists of microscopic heat resistor elements, which convert the electrical energy into 
heat via the Joule effect. The electric pulses thus converted into thermal signals manifest themselves as heat transferred 
to the surface of the thermographic material, e.g. paper, wherein the chemical reaction resulting in colour development 
takes place. This principle is described in "Handbook of Imaging Materials"' (edited by Arthur S. Diamond - Diamond 
Research Corporation - Ventura, California, printed by Marcel Dekker, Inc. 270 Madison Avenue, New York, ed. 1 991 , 

25 p . 498-499). 

[0005] A particular interesting direct thermal imaging element uses an organic silver salt in combination with a re- 
ducing agent. An image can be obtained with such a material because under influence of heat the silver salt is developed 
to metallic silver. 

[0006] Referring to figure 1 , there is shown a global principle schema of a thermal printing apparatus 1 0 that can be 

30 used in accordance with the present invention (known from e.g. EP 0 724 964, in the name of Agfa-Gevaert). This 
apparatus is capable of printing lines of pixels (or picture elements) on a thermographic recording material m, com- 
prising thermal imaging elements or (shortly) imaging elements, often symbolised by the letters le. As an imaging 
element le is part of a thermographic recording material m, both are indicated in the present specification by a common 
reference number 5. The thermographic recording material m comprises on a support a thermosensitive layer, and 

35 generally is in the form of a sheet. The imaging element 5 is mounted on a rotatable platen or drum 6, driven by a drive 
mechanism (not shown) which continuously advances (see arrow Y representing a so-called slow-scan direction) the 
drum 6 and the imaging element 5 past a stationary thermal print head 20. This head 20 presses the imaging element 
5 against the drum 6 and receives the output of the driver circuits (not shown in Fig. 1 for the sake of greater clarity). 
The thermal print head 20 normally includes a plurality of heating elements equal in number to the number of pixels in 

40 the image data present in a line memory. The image wise heating of the heating element is performed on a line by line 
basis (along a so-called fast-scan direction X which generally is perpendicular to the slow-scan direction Y), the "line" 
may be horizontal or vertical depending on the configuration of the printer, with the heating resistors geometrically 
juxtaposed each along another and with gradual construction of the output density. Each of these resistors is capable 
of being energised by heating pulses, the energy of which is controlled in accordance with the required density of the 

45 corresponding picture element. As the image input data have a higher value, the output energy increases and so the 
optical density of the hardcopy image 7 on the imaging element 5. On the contrary, lower density image data cause 
the heating energy to be decreased, giving a lighter picture 7. 

[0007] In the present invention, the activation of the heating elements is preferably executed pulse wise and prefer- 
ably by digital electronics. Some steps up to activation of said heating elements are illustrated in Figs 1 and 4. First, 

so input image data 1 6 are applied to a processing unit 1 8. After processing and parallel to serial conversion (not shown) 
of the digital image signals, a stream of serial data of bits is shifted (via serial input line 21 ) into a shift register 25, thus 
representing the next line of data that is to be printed. Thereafter, under control of a latch enabling line 23, these bits 
are supplied in parallel to the associated inputs of a latch register 26. Once the bits of data from the shift register 25 
are stored in the latch register 26, another line of bits can be sequentially clocked (see ref. nr. 22) into said shift register 

55 25. A strobe signal 24 controls AND-gates 27 and feeds the data from latching register 26 to drivers 28, which are 
connected to heating elements 29. These drivers 28 (e.g. transistors) are selectively turned on by a control signal in 
order to let a current flow through their associated heating elements 29. 

[0008] The recording head 20 is controlled so as to produce in each pixel the density value corresponding with the 
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processed digital image signal value. In this way a thermal hard-copy 7 of the electrical image data is recorded. By 
varying the heat applied by each heating element to the carrier, a variable density image pixel is formed. 
[0009] A control algorithm must determine for every heating element the amount of energy which must be dissipated. 
In practice, the controller algorithm must deal with a variety of real-world problems: 

5 

o Changing characteristics of the film media give different pixel sizes for the same nib (or heating element) energy, 
e.g. some examples: 

a different humidity in the emulsion layer, making its thermal capacity different, 
10 - a different chemical composition of the image forming components. 

* Environmental characteristics like temperature and humidity may change: 

a temperature rise of the environment must be taken into account as the image forming temperature will not 
15 rise and is given by the chemical composition of the emulsion layer, 

humidity again changes the thermal capacity of the emulsion, producing different temperature rises when 
applying the same amount of energy. 

o The thermal process itself produces an excessive amount of heat which is not absorbed by the image forming 
20 media. This excessive heat is absorbed by a heat sink, but nevertheless, gives rise to temperature gradients 

internally in the head, giving offset temperatures in the nibs and between the several nibs. E.g. when the image 
forming process must have an accuracy of 1 °C in the image forming media, an increased offset temperature of 
5 °C in the heat generating element must be taken into account when calculating the power to be applied to that 
element. 

25 © The heat generating elements are in the ideal case fully thermally isolated from each other. In practice, this is 
never the case and cross-talk exists between the several nibs. This cross-talk can be localised on several levels: 

heat transfer between the several nibs in the thermal head structure itself, 
heat transfer in the emulsion and film layer itself, 
30 - pixels are not printed one aside the other, but partly do overlap on the print media, mechanically mixing heat 

from one pixel with the other. 

o The electrical excitation of the nibs is mostly not on an isolated base. This means that not every nib resistor has 
its own electrical voltage supply which can be driven independent of all the other nibs. In general, some drive 

35 signals for driving the nibs are common to each other, this with the purpose of having reduced wiring and drive 

signals. In general, all nibs can be only switched on or off in the same time-frame. Producing different weighted 
excitations can only be achieved by dividing the excitation interval in several smaller intervals where for every 
interval, it can be decided if the individual nib has to be switched on or off. This process of "slicing" has its influence 
on the thermal image forming process. For example, a pattern excitation with the weights (1 28,0,0,0,0,0,0,0) differs 

40 mathematically only 1 point from a pattern excitation with the weights (0,64,32,16,8,4,2,1), but the pixelsize often 

will be much more different than just 1 point. It has been perceived that in some thermal heads, even a 'zero- 
excitation' or a 'no excitation' interval produces some heat in the nib as well. The controller must take this effect 
into account. 

45 [0010] An empirical way of trying to solve the mentioned problem could comprise the steps of making a printout of 
all the available slices, measuring the density or pixel size on the pixel output, and deriving a relation between pixel 
output and the slices used. By simply using a conversion table one could build a continuous monotone and maybe 
even a linear relation between my table index and the pixel output. 
However, such method is not feasible, because of several reasons, among which: 

50 

° Only the large excitation times (or slices) will give a pixel output. 

o Also the smaller excitation times are important, as they are still used for compensations, even though no pixel 
output can be detected on the thermographic material. 

• Density measurement or pixel size measurement is not always error free, making the evaluation of the results 
55 difficult and asking more for a statistical evaluation of the results. 

° When making a printout, a temperature rise will occur in the printing device, which can jeopardise the whole meas- 
urement. 
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[0011] Although it is known to prepare both black-and-white and coloured half-tone images by the use of a thermal 
printing head, a need for an improved recording method still exists. 

OBJECTS OF THE INVENTION 

5 

[0012] It is an object of the present invention to provide an improved method for recording an image on a thermal 
imaging element by means of a thermal head having energisable heating elements. 

[0013] Other objects and advantages of the present invention will become clear from the detailed description and 
the drawings. 

10 

SUMMARY OF THE INVENTION 

[0014] The above mentioned object is realised by a method for generating an image on a heat mode imaging element 
having the specific features defined in the independent claims. Specific features for preferred embodiments of the 
15 invention are disclosed in the dependent claims. 

[0015] Further advantages and embodiments of the present invention will become apparent from the following de- 
scription and drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

[0016] The invention will be described hereinafter with reference to the accompanying drawings (not necessarily to 

scale), which are not intended to restrict the scope of the present invention. 

Herein, 

25 Figure 1 shows some basic functions of a direct thermal printer; 

Figure 2 shows a function f1(x) which is continuous in all points, but which is not continuous in point a; 

Figure 3 shows a function f2(x) witch is continuous over the full region but not monotone in the region b-c; 

Figure 4 shows a control circuitry in a thermal print head comprising resistive heating elements; 

Figure 5 shows a side view of a thick film head in contact with a thermographic material; 
30 Figure 6 shows an enlarged cross section of a thick film head and the different material layers; 

Figure 7 is a graph illustrating a parallel to serial conversion; 

Figure 8 shows a graphical representation of the function u(t,t start) At); 

Figure 9 shows the power applied to a heating element for a sequence with maximum weight; 

Figure 10 shows the heating element power supply for a sequence {1 ,0,1 0}; 

35 Figure 1 1 is a chart illustrating for all heating elements the activating heating pulses with an exemplary duty-cycle 

and with an exemplary skipping, known from prior art; 

Figure 12 shows different temperature rises in a heating element for two equal energy excitations; 
Figure 13 shows a macro density of a line pattern as a function of an increased electrical excitation for binary 
weighted time slices, as suitable in a method according to the present invention; 
40 Figure 1 4 shows a 1 D model for a thermal head structure, in a line running from the heatsink, through the nib and 

the thermographic material up to a platen , as suitable in the present invention; 

Figure 15 shows a maximum temperature in a thermographic material as a function of an increased electrical 
excitation weight, as suitable in a method according to the present invention; 

Figure 16 shows the measured pixel size from a macro density measurement for a binary slicer with weights in 
45 the range from 125 to 255, as suitable in a method according to the present invention; 

Figure 1 7 shows the calculated T max function for binary slices with weights in the range from 1 25 to 255, as suitable 
in a method according to the present invention; 

Figure 18 shows calculated T max values for a slicer with weight vector {1,2,4,8,16,32,32,32} and several points 
marked for inclusion to the LUT2 table, as applicable in a method according to the present invention; 
50 Figure 19 shows the theoretically calculated maximum temperature Tmax in a thermographic material m for a 

continuous and monotone set of sequences according to the present invention; 

Figure 20 shows the measured pixel sizes on a thermographic material for a continuous and monotone set of 
sequences according to the present invention. 

55 DETAILED DESCRIPTION OF THE INVENTION 

[0017] The description given hereinafter mainly comprises seven sections, namely (i) terms and definitions used in 
the present application, (ii) general description of a preferred embodiment of a method for thermal printing according 
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to the present invention, (iii) preferred embodiment of a method for thermal printing comprising an experimental building 
of a look-up table, (iv) further preferred embodiment of a method for thermal printing according to the present invention, 
(v) illustration of the invention with practical examples, (vi) further preferred embodiments according to the present 
invention, and (vii) further applicability of a method according to the present invention. 
5 [0018] While the present invention will hereinafter be described in connection with preferred embodiments thereof, 
it will be understood that it is not intended to limit the invention to those embodiments. 

(i) Explanation of terms used in the present description 

10 [001 9] For the sake of clarity, the meaning of some specific terms applying to the specification and to the claims are 
explained before use. Most of these terms are explained in this section, some more specific terms [e.g. af or u(t,t start , 
At) or p H ^f)\ 8k ] will be explained in the appropriate paragraphs (e.g. while explaining equation 5). 
[0020] An "original" is any hardcopy o r softcopy containing information as an image in the form of variations in optical 
density, transmission, or opacity. Each original is composed of a number of picture elements, so-called "pixels". Further, 

is in the present application, the terms pixel and dot are regarded as equivalent. 

[0021 ] Furthermore, according to the present invention, the terms pixel and dot may relate to an input image (known 
as original) as well as to an output image (in softcopy or in hardcopy, e.g. known as print). 

[0022] The term "thermographic material" (being a thermographic recording material, hereinafter indicated by symbol 
m) comprises both a thermosensitive imaging material and a photothermographic imaging material (being a photosen- 
20 sitive thermally developable photographic material). 

[0023] For the purposes of the present specification, a "thermographic imaging element" le is a part of a thermo- 
graphic material m (both being indicated by ref. nr. 3). 
Hence, symbolically: m 3 le. 

[0024] By analogy, a thermographic imaging element le, comprises both a (direct or indirect) thermal imaging element 
25 and a photothermographic imaging element, in the present application the term thermographic imaging element le will 
mostly be shortened to the term imaging element. 

[0025] By the term "heating material" (hereinafter indicated by symbol hm) is meant a layer of material which is 
electrically conductive so that heat is generated when it is activated by an electrical power supply. 
[0026] In the present specification, a heating element H n is a part of the heating material hm . 
30 Hence, symbolically: hm 3 H n . 

[0027] A "heating element H n " (as e.g. H 1 , H 2 , H 3 ...; sometimes also indicated as "nib") being a part of the heating 
material hm is conventionally a rectangular or square portion defined by the geometry of suitable electrodes. 
[0028] A "platen" comprises any means for firmly pushing a thermographic material against a heating material e.g. 
a drum or a roller). 

35 [0029] According to the present specification, a heating element is also part of a "thermal printing system", which 
system further comprises a power supply, a data capture unit, a processor, a switching matrix, leads, etc. 
[0030] The index "n" is used as an superscript with regard to nib numbers, n=0,1 ,...,N nibs -1 with N nibs the total number 
of nibs on the thermal head. 

[0031] The term "controllability" of a thermal printing system denotes the ability to precisely control the output of a 
40 pixel, independent from the position of the pixel and the presence of pixel neighbours. 

[0032] A "heat diffusion process" is a process of transfer of thermal energy (by diffusion) in solid materials. 

A "heat diffusion partial differential equation PDE" is a partial differential equation describing a heat diffusion process 

in a solid material. 

A "specific heat production q n " is a volumetric specific thermal power generation in the confined bulk of the thermo- 
45 graphic material [W/m 3 ]. When reference is made to a sequence S k defining the power distribution, q k n will be used. 
A "specific mass p" is a physical property of a material and means mass per volumetric unit [kg/m 3 ]. 
A specific heat c means a coefficient c describing a thermal energy per unit of mass and per unit of temperature in a 
solid material at a temperature T [J/kg.K] 

[0033] A "thermal conductivity V is a coefficient describing the ability of a solid material to conduct heat, as defined 
so by Fourrier's law 



55 X is expressed e.g. in [W/(m.K)]. 

[0034] A "transient temperature history" Th n = f(x,y,z,t), or short Th n , denotes the temperature history in the thermo- 
graphic material in every point x,y,z at pixel n and at the instance of time t (t being relative to the start time of the heating 
element activation). Th k n or Th k refers to a sequence S k or power distribution P k used as a template for driving the 
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heating element. 

[0035] Dh k will be denoted as a calculated quantity from a theoretical background and being representative for the 
considered pixel output D k for a given temperature history Th k in the image forming material. 

[0036] GS denotes a set of Dh k values found by calculating the pixel output for each given excitation S k or P k> GS= 
5 {{S 1 ,Dh 1 ) l (S 2 ,Dh 2 ), (S 3l Dh 3 ),...} orGS^P^Dh,), (P 2 ,Dh 2 ), (P 3 ,Dh 3 ),...}. 

[0037] GS' and GS" each denote an ordered subset from GS that imposes a mathematical relationship on the Dh k 
values. 

[0038] An "activation pulse" is an energy pulse supplied to a heating element, described by a certain energy given 
during a defined time interval ts. 
10 The elementary time interval ts during which a strobe signal is active is often called a "time slice". 

The term 'lime slice s, of activation pulses" explicitly indicates that during a time slice, and hence during a same strobe 
signal, the individual heating elements may be individually and independently activated or non activated by correspond- 
ing activation pulses. 

[0039] The "index i" is used to denote a slice number, i=0,1 N-1 , with N being the total number of slices. 

15 [0040] A time slice ts may be written as 

Equation 1 : ts = tScale * w 

20 wherein w is the weight of the slice and wherein tScale is a scaling factor which is the same for all slices. 

The term "weight w of a slice" indicates a relative relationship among the time slices used. The smallest slice weight 
is mostly taken as being 1 . 

The term "binary weights" refers to a set of weights wherein the increase in magnitude is a power of 2, symbolically Wj 
=2i, wherein je{0,1 ,2,...,i I ...,N-1} and each number picked only once. 
25 [0041] A "sequence S k " means a combination of time slices with a particular order and each with a particular weight 
or pulse length ts that can be different for each time slice (as is illustrated e.g. in Figs. 9 and 10). 
The "index k" is used to denote sequences of slices, k=0,1 ,.,.,K-1 , with K the total number of possible sequences or 
the total number of possible combinations of slices. 

[0042] The "weight W k " of a sequence S k is the sum of the weights of the active slices in that sequence. 
30 The total excitation time t k n for a nib n due to a sequence Sk equals y 1 =t Sca | e * W k . 

[0043] The term "sequence S k of time slices Sj having (different) binary weights Wj" comprises a collection of time 

slices using a binary weighting scheme, wherein tSj is different for each slice i and tSj=tScale*2J, j being associated with 

slice number i but not necessarily applied in numerical order or je{0,1 ,2 1 N-1} and each number picked only once. 

[0044] The term power distribution P k or P k (t) is denoted for a time dependent power cycle applied to a heating 
35 element during a line time and is more general than a sequence of slices, mostly supplying the heating element with 

a constant voltage. P k is a more general formulation of a heating element excitation, the index k referring to a number 

identifying the power distribution, k=0,1 ,...,K-1 . 

[0045] In the present application, a "pixel output D" or shortly an "output D" comprises a quantification of a pixel 
printed on a recording material, said quantification possibly relating to characteristics as density, size, etc. The pixel 
40 output of nib n is denoted D n and the pixel output for a particular sequence S k or power distribution P k is denoted D k . 
.[0046] A "continuous output" means an output showing a continuous output for a continuous rising of the input to 
the system. 

The term "continuous" has to be interpreted in a mathematical sense. Therefor, for a function f(x), rt applies for every 
co-ordinate that 

45 

Equation 2 : f(a) = lim /(jc) = lim f(x) . 

jr— *a* x—*a~ 

50 

[0047] A "monotone output" means an output showing a monotone output for a monotone rising of the input to the 
system. [In the context of this patent, the input preferably is the index of a table LUT, representing in fact the input to 
a virtual printing device, making abstraction of the time slice details.] 
The term "monotonous" can mathematically be written as: 

55 

Equation 3: Vx,y e x > y => /(*) > f(y) . 
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The acronym M D k " represents a set of outputs with no particular ordering or relation between the elements with regard 
to the index k. 

The acronym "D kcm " represents a set of ordered pixel outputs with an order being defined by a mathematical relationship 
between the elements, in this case imposing continuity and monotony with regard to the index k. 
[0046] The function F cm (k) denotes a function being continuous and monotone with regard to the input k. It is an 
application defined curve and can take whatever mathematical form, as long as it is continuous and monotone. 
[0049] The acronym M D kcm |" represents a set of ordered pixel outputs with an order being defined by a mathematical 
relationship between the elements, in this case imposing continuity, monotony and linearity with regard to the index k, 
or D kcmi = ( r i * k) + r 2 , wherein ^ and r 2 are real numbers. 

[0050] The function F cm ,(k) is a function being continuous, monotone and linear with regard to the index k. It is an 

application defined curve in the form F cm ,(k) = (^* k) + r 2 , wherein r, and r 2 are real numbers. 

[0051] An output D k is related to a sequence S k . Mathematically a formal notation can be used, D k = F(S k ), implying 

that a relation exists between a sequence S k of slices of activation pulses and a pixel output D k . 

[0052] A lookup table LUT1 can be defined that relates a certain index number k to a certain pixel output D k or Dh k 

in a first entry and in a second entry the sequence S k or the power distribution P k to be used. Based on LUT1 , a relation 

can be established between the index k of the LUT and a pixel output D k (or Dh k ), further related to a sequence S k or 

power distribution P k , to be used in the printing process. No particular mathematical relationship is necessary between 

the index k of the table and the numerical value of D k or Dh k . It is a numerical representation of the set GS. 



Index 


Wanted pixel output 


Excitation sequence 


1 


D 1 or Dr^ 


S 1 or P 1 


2 


D 2 or Dh 2 


S 2 or P 2 


k 


D k or Dh k 


S k or P k 



[0053] A lookup table LUT2 can be defined that relates a wanted pixel output D kcm or Dh kcm to a sequence S k or 
power distribution P k to use. Furthermore, with regard to the index k of the table, the D kcm or Dh kcm values will be a 
continuous and monotone set of numbers. LUT2 is the numerical presentation of the set GS'. 

[0054] A lookup table LUT3 can be defined that relates a wanted pixel output D kcml or Dh kcm , to a sequence S k , or 
a predefined time dependent power distribution P k to use. Furthermore, with regard to the index k of the table, the 
D kcm | or Dh kcm] values will have a continuous, monotone and linear relationship with the index k. LUT3 is the numerical 
presentation of GS". 

[0055] A vector of slice weights will be denoted by {w 0 , w 1f w 2 w N } with w 0 the weight used for the first slice, w 1 

for the second slice etc.... 

[0056] As mentioned before, some more specific terms [e.g. a, k or u^^^At) or Pud^s*] w '" be explained in the 
appropriate paragraphs. 

[0057] It is known, and put to intensive commercial use (e.g. Drystar ™, of Agfa-Gevaert ), to prepare both black- 
and-white and coloured half-tone images by the use of a thermal printing head, a heat-sensitive material (in case of 
so-called one-sheet thermal printing) or a combination of a heat-sensitive donor material and a receiving (or acceptor) 
material (in case of so-called two-sheet thermal printing), and a transport device which moves the receiving material 
or the donor-acceptor combination relative to the thermal printing head. In the next paragraphs, a working method 
according to the present invention will be explained in full depth. 



(ii) General description of a preferred embodiment of a method for thermal printing according to the present invention. 



50 



55 



[0058] According to the present invention, a method for printing an image using a thermal printing system comprises 
a thermal printer 10 having a thermal head TH (ref. nr. 20) incorporating a plurality of energisable heating elements 
H n (29) and a thermographic material m (5), said method comprising the steps of 

supplying image data Id (see ref. nr. 16) corresponding to a plurality of output values D k to a processing unit 18 
of said thermal printer, 

reading a first table LUT1 comprising first entries for desired output values D k and second entries for sequences 
S k of time slices s, of activation pulses, 

transformation of said image data Id into corresponding sequences S k of time slices of activation pulses, 
providing said time slices of activation pulses to said heating elements H n of said thermal head TH, 
printing said image by transporting (see ref. nr. 8 and arrow Y) said thermographic material past and adjacent to 
said thermal head and by activating (see ref. nr. 30) said heating elements H n of said thermal head, 
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wherein said transformation of said image data Id comprises a step of establishing a relation between output values 
D kcm which are both continuous and monotone, and corresponding sequences S k of time slices s, of activation pulses. 
[0059] In the following paragraphs, each of these method-steps will be described in full detail. 
The step of supplying image data l d corresponding to a plurality of output values D k to a processing unit of said thermal 
printer has already been indicated before in relation to the description of Figure 1, in particular w.r.t. input data block 
16. Here, forsake of good understanding the present invention, some further informations are explained with regard 
to Figs. 7, 8, 9, 10 and 11. 

[0060] First, we refer to figure 7 which is a graph illustrating a parallel to serial conversion. Knowing that the thermal 
head normally includes a plurality of heating elements equal in number to the number of pixels in the data present in 
the line memory and that each of the heating elements is capable of being energised by heating pulses, the number 
of which is controlled in accordance with the required density of the corresponding picture element, Fig. 7 illustrates 
the conversion of a ten-head-row subjected to image data of bytes consisting of two bits, and thus representing max- 
imally four densities (image data '00' gives a printing sequence {000}, '01* gives {100}, '10' gives {110} and '11' gives 
{111 }). Herein, following referral numbers are used: 40 is a schematic for a parallel to serial conversion, 42 is schematic 
for serial outputs, and 44 is a schematic for parallel inputs. It follows that the thermal head applied with a recording 
pulse causes current to flow through the heating elements for corresponding "ones" (cf. input data indicative of "black 
picture elements") of. 

[0061] Fig. 8 is the basic function used to describe the strobe signal 50 applied to the electronics (by means of strobe 
Iine24), being active during a limited time At starting from t = t^^. During a line time, N strobe pulses will be successively 
applied to the system, where At can be different for each strobe pulse (as is illustrated e.g. in Figs. 9 and 10). We will 
denote the word slice to each sequence where the strobe pulse is active. During a single time slice, the power applied 
to a heating element preferably will be fixed. 

[0062] Fig. 9 shows the power applied to a heating element for a sequence with maximum weight, and describes a 
sequence of N successive strobe pulses applied to the heating element during a line time. In this case, during each 
slice, the heating element will get a power P on . The power sequence applied to the heating element can be represented 
by a sequence of numbers {aj that represent the power state applied to the heating element during each slice, '1 1 being 
P on and '0' being P off applied to the heating element. For the depicted example, all a r values are 'V or {a,}={1,1,1 
1 }. When N slices have been given, the power to the heating elements is switched off. For sake of clarity, it is indicated 
that ref . nr. 52 is a time slice having a weight iv 0 and being in an "on -state" , 54 is a time slice having a weight w, and 
being in an "on-state" , and 56 is a time slice having a weight vv NnJbr ^ and being in an "on-state". 
[0063] Fig. 1 0 shows the heating element power supply for a sequence {1 ,0,1 ,...,0}, and hence describes a sequence 
of N slices with only P on applied in the first and third slice and P off during all the other slices. This can be represented 
by { a iM1 , 0, 1 , 0, 0,..., 0}. Ref. nr. 58 represents a time slice having a weight w A and being in an 11 off -state". 
[0064] Figure 11 is a 3D- chart 60 illustrating strobe pulses with an exemplary but fixed duty-cycle and with an ex- 
emplary skipping) for all heating elements (up to nib nr n), as known from prior art. In Figure 11 , skipped pulses are 
indicated by dotted lines and following referral numbers are used: 61 represents a first time slice for all nibs, 62 is 
second time slice, 63 is a last time slice, 64 is a first activation pulse for a first nib 0 , 65 is a first activation pulse for a 
last nib N nibs -1 , 66 is the last activation pulse for first nib 0 , 67 is the last activation pulse (and being skipped) for last 
niD N nibs" 1 . 68 is an active activation pulse (T) applied to the last-but-one nib N njbs -2 during the third slice, 69 is a 
skipped activation pulse or powerless activation pulse ('0') applied during the third time slice. 

[0065] The method step of reading a first table LUT1 comprising first entries for desired output values D k and second 
entries for sequences S k of time slices s, of activation pulses, will be disclosed hereinafter. 

• The technique of reading a table, e.g. a look-up-table, comprising first entries and second entries, is generally well 
known to people skilled in the art and hence, does not need further explanation. In this particular table LUT1 , the 
first entries relate to desired output values D k , and the second entries relate to sequences S k of time slices s, of 
activation pulses. The terms 'sequences S k of time slices s, of activation pulses' have already been indicated in 
foregoing section (i) related to "terms and definitions" but still will be explained more thoroughly in the next para- 
graphs. 

• In order to be as clear as possible, first attention is focused on the electric driving of the individual heating elements 
which is based on the schematic shown in Figure 4. Essential is the presence of a STROBE signal 24 which is 
common for a whole series of heating elements 29 and which allows a common on- or off-switching of all heating 
elements. 

• In order to be able to give the heating elements 29 a different electrical energy (e.g. because of compensation 
reasons), and taking into account that several electrical signals are common to all heating elements, a so-called 
"time slicing" or shortly "slicing" is performed. Said slicing comprises a dividing of the total excitation period in 
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shorter individual time periods or "slices" and selectively activating heating elements during each time slice (by 
means of AND-gates ref. nr. 27). 

© More in particular, the total excitation time is divided into small time slices, and the total activation time of a heating 
s element is then produced by turning the heating element on or off in some of those slice intervals. Pure arithmet- 

ically, a total excitation time equals the summation of the on-times of the different slices and is in the end nearly 
equal to the desired excitation time for the heating element. So, slices with weights of e.g. 1 ,2,4,8,16,32,64,128 
allow to make any random weight between 0 to 255. 

The step of transforming said image data tointo corresponding sequences S k of time slices of activation pulses, 
10 may be explained by the following. 

o As mentioned before, each of the heating elements is capable of being energised by activation pulses, the energy 
of which is controlled in accordance with the required output value (e.g. density, or pixel size) of the corresponding 
pixel. 

15 

° When image data have a relatively high value, the corresponding output energy is also relatively high, and so is 
e.g. the optical density of the hardcopy Image 7 on the thermographic material 5. On the contrary, lower density 
image data cause the heating energy to be decreased, giving a lighter picture 7. 

20 © The output energy E nftJ dissipated in a heating element switched on during a time t^ is representative for the final 
temperature in the thermographic material m and equals in its most simply form for a resistive heating element 29: 



v 2 

Equation 4: E nib = D *™ t a 



25 ' " on 



[0066] Herein, the term V,^ represents the voltage applied to the heating element during an on-time, 
R heatereiement tne equivalent macroscopic resistance measured between the terminals supplying the heating element, 
t on represents the sum of the time slices where the heating element was switched on. 
30 [0067] The step of providing said time s/ices of activation pulses to said heating elements H n of said thermal head 
TH needs no further explanation, as it is well known by people skilled in the art. 

[0068] For a same reason, also the step of printing said image by transporting said thermographic material past and 
adjacent to said thermal head (comprising so-called fast scan movement X and a slow scan movement Y), and by 
activating said heating elements H n of said thermal head, needs no further explanation. 
35 [0069] Now an important characteric will be explained, in particular the characteristic wherein said transformation of 
said image data id comprises a step of establishing a reiation between output values D kcm which are both continuous 
and monotone, and corresponding sequences S k of time slices s } of activation pulses. 

[0070] According to the present invention, we have found that the output D, especially D kcm , is dependent from the 
time history of activation, in particular dependent from the history of time slices s ( of activation pulses, and more in 

^o particular dependent from the sequences S k of time slices s, of activation pulses. 

[0071] In simple language, it is said that time history is an important parameter in the present invention, as it is dealing 
with the transfer of thermal energy from inside the heating element towards the thermographic material. Temperature 
behaves as a non-linear function of time. For example, if 1 mJ heat is produced in the heating element, followed by a 
300 \xs wait period and then again a 1 mJ heat production, this will have a different thermal effect on the thermographic 

45 material as giving a single 2 mJ single heat pulse. In the former case, most of the heat of the first 1 mJ pulse will have 
vanished before the influence of the second heat pulse can be noticed. A calculated example is shown in Figure 12. 
Although equal energy has been given to the heating element, the sequence of giving this energy is different and this 
gives a different history of the temperature. When giving the energy in a single pulse (ref. nr. 72), a much higher 
temperature is reached, giving in practice a more pronounced pixel output than the split energy scheme (ref. nrs. 74, 76). 

so [0072] In a more mathematical wording, it can be said that a general expression of the time dependent power applied 
to the heating element using a time slice system: 



Equation 5 : P„e(OL n(f.C.w, t 5eale )[a* P on + {\-a-)P off ) 
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[0073] For sake of easy reading, the abovementioned nomenclature (see section i) further comprises: 

Ptidfysk^e time dependent power applied to a heating element when using a sequence S k , 
P on is the power dissipated in a heating element when it is active driven, 
s P off the power dissipated in the heating element when it is switched off during an active write period, i.e. an active 

STROBE pulse, 

u(t,t start ,At) a function that is T starting from the time during a time At (see Figure 8), 
f the active starting time of the i'th slice, 
w[1s the weight of a slice s j( 
10 N is the number of slices, 

i is a summation index (over all time slices) 

aj k is a sequence coefficient having only the values of 0 or 1 and defines whether the heating element is on (1 ) or 
off (0) during the slice i, the index k referring to the index of the sequence S k , 

the aj k values, together with the order in which they appear will be called a sequence S k and are therefore closely 

is linked to the defined slices, {a^a^ a N-1 k }, 

a slice is a time span ts, having an heating element activation time defined by tSj= w, * t Sca , e . 

[0074] For an image forming process using a slicing method, the total excitation time of a heating element is distrib- 
uted over a large time period. During some slices, a heating element will have a { equal to zero, making it cool down 
20 during that period. For some thermal heads, it is also possible that an a, value of zero, still gives some limited power 
excitation in the heating element itself during an active slice period. 

[0075] Now we have found that the output D k , especially D kcm is not only dependent from the energy of activation, 
but in particular dependent from the time slices s i of activation pulses, and more in particular dependent from the 
sequences S k of time slices s, of activation pulses. So it may be that the application of a energy E., to a particular 
25 heating element H n of a sequence S 1 of time slices s M of activation pulses results in a different output as compared to 
the application of an equal energy E 2 = Ej to the same heating element H n of another sequence S 2 of time slices s i2 
of activation pulses. 

[0076] During each slice period i, the heating elements will be STROBE active during a time ts { = w s * t Scale and 
depending on the value of aj, the heating element will be put on or off. In this manner, random values of t n can be 
30 realised by choosing a smart vector of time weights, a sufficient number of slices N and a correct time scaling factor 
t Scate . It is important to understand that a same excitation time t n can be made in different ways, as it is always possible 
to change the order of the successive slices, giving in total N! combinations. 

[0077] If we call t n the time that the heating element n must be excited by electrical power, then we have for printing 
a line, a whole vector {t n } with different t n values for the several heating elements, representing the image data. A well 
35 known way of slice weights is e.g. a binary series, having the weights {1 ,2,4,8,1 6,32,64,...}. Depending on the maximum 
value of t n and the sum of all the weights, the correct scalingf actor t^^ can be chosen. 

[0078] When using a sequence of slices, it is often assumed that the pixel output shows a continuous and monotone 
"pixel output" (e.g. visual density or pixel size) in accordance with the used excitation weight W k , being representative 
for the integrated power over all slices of the sequence, and in case of a constant voltage excitation, linear with the 
40 total on-time t k n = W k * tg^. In practice, the time dependent character of the power distribution applied to the heating 
element, has a substantial influence on the pixel output produced. When using a series of sequences with slices where 
the summed excitation time (t k n ) increases continuously, often a non-monotone pixel appearance (e.g. density of pixel 
size) will result. 

[0079] In practice, a method according to the present invention calculates the temperature response for all the avail- 
45 able sequences. Then, based on the calculated pixel output, one determines slice sequences out of the pool of available 
sequences and builds a set of sequences giving a continuous, monotone -and maybe even linear- pixel output when 
applied to the head. This data are put in a table, having as input a desired pixel size or density, and as output a sequence 
describing a slice based power distribution to be applied to the heating element. 

50 (iii) Preferred embodiment of a method for thermal printing comprising an experimental building of a look-up table. 

[0080] First, it is important to know that a small heating element on top of a substrate produces the heat for interacting 
with a thermographic material. Only a small part of the heat will flow towards the thermographic material, the rest are 
losses which will flow through the substrate to a heatsink structure, kept at a low temperature by means of cooling 
ss techniques. 

[0081] Second, it is also important to know that different heating elements on a same thermal head generally get 
different amounts of energy during a same print cycle, because of several possible reasons which are not explained 
here because of conciseness. In a practical experiment, a periodic line pattern has been printed on a thermographic 
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material with an increasing value of the energy sent to the heating elements according to the above mentioned total 
on-time t" of the heating element n. The weights used where {1 ,2,4,8,16,32,64}. By means of a macro densitometer, 
the density of the periodic line pattern has been recorded, as this is sufficiently representative for the line thickness. 
[0082] Here, reference is made to figure 1 3 showing a macro density of a tine pattern as a function of an increased 

5 electrical excitation time for binary weighted time slices. Herein, 80 represents a density curve relating to sequence 
weights W k , 82 a non linearity in density at a sequence weight of 64, ref. nr. 84 a jump in density at a sequence weight 
of 128, ref. nr. 86 a jump in density at a sequence weight of 196, and ref. nr. 88 represents a saturation in density. 
[0083] Figure 1 3 shows clearly that for some weight transitions (cf. x-axis giving W k of sequence), important changes 
in the pixel output exist. An important jump in the weight distribution W k over time occurs at 1 28, 1 96, 64. 

10 [0084] These jumps can be removed by using a lookup table and using the concept of sequences S k . The input of 
such table is a virtual excitation and is representative for the output result from the controller. The output of the table 
is a sequence S k that actually has to be used for driving the heating element. The values in the table can be picked 
on an experimental base by making pixel output examples and then choosing those sequences S k so that a continuous 
and monotone result is obtained. As an example, taking Figure 13, the first 1 28 values of our table could be the values 

15 from 1 to 128, as these weights give an acceptable continuous and monotone output. Entry number 129 should not 
have weight 129, as this gives a major drop in our pixel output (e.g. density drops from 0.4 to 0.33). The first weight 
having a density a little higher than 0.4 is about weight 155. So, this becomes the value in the lookup table at entry 
number 129. For this example, it is possible that one would not be able to construct a table with 255 different "grey- 
values", but somewhat less. 

20 [0085] Nevertheless, when applying the above method for a practical thermal head on a certain thermographic ma- 
terial, some problems do arise making the construction of the lookup table very difficult. 

In general, the experimental way to make a lookup table is possible but is expensive because of the many experiments 
that have to be done under well controlled conditions. 

25 (iv) Detailed description of a further preferred embodiment of a method for thermal printing according to the present 
invention. 

[0086] For sake of good understanding, first some general introduction is given relative to this preferred embodiment. 
[0087] Concerning the image forming process, it is assumed that a model or calculation scheme is known, giving e. 
30 g. a relation between the temperature history in a point of the thermographic material and its corresponding pixel output 
expressed e.g. in a density variation or corresponding pixel size. This model includes the time dependency of the heat 
applied to the thermographic material. 

In case the correct details of this formulation are not known exactly, a simple linear relation may be taken, e.g. pixel 
size varying linearly with peak temperature in the thermographic material. The exact values of this relation are not 
35 important, only the relative behaviour is important, as the most important thing is to get a continuous monotone relation 
between excitation weight and pixel output. 

[0088] More in particular, the invention consists of making a finite-element-model FEM or a finite-difference-model 
FD calculation scheme for the head, including heatsink at a constant temperature, isolation layers, heating element, 
protective layer, thermographic material and the material pressing the film against the thermal nib line (mostly rubber). 

40 Using the known way in which the heating element is excited for the given sequence of slices, a transient simulation 
is made of the thermal partial differential equations PDE governing the thermal system. This simulation gives as output 
the history of the temperature in the thermographic material. Once the temperature is known, the response of the pixel 
output can be found for the given power distribution applied to the heating element, defined by the sequence S k of slices. 
[0089] In a further preferred embodiment of a method for thermal printing according to the present invention, said 

45 step of establishing a relation between output values D kcm which are both continuous and monotone, and corresponding 
sequences S k consisting of time slices s, of activation pulses, comprises the following substeps: 

building a calculation scheme representative for said thermal printing system, said model taking into account ther- 
mal characteristics of said thermal head TH, thermal characteristics of said thermographic material m, thermal 
50 characteristics of a platen and specific heat productions q k as imposed by sequences S k of said time slices s i 

having different binary weights w ; ; 

calculating transient temperature histories Th k in said thermographic material m at nib n for said sequences S k of 
time slices Sj having different weights w } (e.g. binary weights); 

establishing a relationship between said transient temperature histories Th k and said output values Dh k ; 
55 . storing first entries for continuous and monotone output values Dh kcm and second entries for corresponding se- 
quences S k into a second table LUT2. 

[0090] In the following paragraphs, each of these method-steps will be described in full detail. 
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[0091] The step reading 'building a calculation scheme representative for said thermal printing system, said model 
taking into account thermal characteristics of said thermal head TH, thermal characteristics of said thermographic 
material m, thermal characteristics of a platen and specific heat productions q k n in nib n as imposed by sequences S k 
of said time slices s, having different weights w,, now will be described immediately, 

[0092] A thermal calculation scheme simulates how the heat is transferred from the heating element to a point inside 
the thermal sensitive material. This is a transient simulation, using a precise image of the time dependent heat gener- 
ated inside the heating element, which is fully defined by the used sequence S k . For most thermal heads (e.g. thick 
film or thin film heads), a one-dimensional model 1 D will be sufficiently accurate for making the calculation. As transient 
simulations need to be performed for all different sequences, the smaller the model (e.g. 1D simulations), the more 
beneficial it is with regard to the calculation time. 

All materials need to be defined by their correct parameters, being thermal conductivity X [W/mK], specific heat c [J/ 
m 3 ], and specific mass p [kg/m 3 ]. 

For the 1 D formulation, the boundary conditions will be dirichlet conditions (thus having a fixed boundary temperature) 
ending at one side into the heatsink, being at a fixed temperature during the simulation and at the other side, it will end 
in the structure which presses the thermographic material against the nib structure, being normally also at a fixed 
temperature. 

[0093] This will be illustrated with an example. Consider a thick film head with a resistor lying on a th in layer of glass; 
a cross section of this configuration is shown in Figs. 5-6. This glass 32 is deposited on a ceramic substrate 34 and 
fixed to a heatsink (not shown). Using e.g. a rubber roller 6, the thermographic material is pressed against the heating 
element 

[0094] From this structure, a thermal calculation scheme can be made based on a finite-element-model FEM or on 
a finite-difference-model FD of this structure. A 1D model is indeed a strong approximation of a real 3 dimensional 
structure. But the nature of the heat diffusion process will be the same for the one-dimensional 1 D, two-dimensional 
2D or three-dimensional 3D problem. Our interest goes to the relative behaviour of the temperature reached in the 
thermographic material with regard to a given sequence S k of heat production in the heating element. Of course it is 
beneficial to use a 3D calculation and will of course give the most correct result, but it can be too expensive and too 
time consuming. 

A thermal process as applied to direct thermography is characterised by following general heat equation PDE: 

Equation 6; pc ^=div(X grad(T))+q(x,y,z,f). 

Herein, the function q(jr,y,z,0 gives the specific heat production in every point of the considered space and is exactly 
defined by the sequences S k applied to all heating elements. 

[0095] A 1 D calculation scheme is made, by taking a line intersection through the nib perpendicular on the heatsink. 
In the model, the cross section is normalised to 1m 2 and heat production in the nib is recalculated with regards to this 
normalised cross section . All dimensions in the 1 D-direction must be taken correctly from the real structure as illustrated 
in Figs. 5, 6 and 14. 

[0096] The boundaries of our calculation scheme are simple dirichlet boundaries, assuming e.g. 0°C on these bound- 
aries. Physically, this is the heatsink or the inside structure of the pressing rollers, both with large time constants, 
allowing a constant temperature in these materials during the simulation. 

[0097] For each material found in the 1 D cross-section, the correct material properties must be inserted in the cal- 
culation scheme 90 (cf. p, X and c; defined hereabove in the section "Explanation of terms used in the present inven- 
tion"). In the example of Fig. 14, we have the index s for a substrate 34 (e.g. a ceramic substrate), the index b for a 
bonding layer 32, the index n for the heating element or nib 29, the index p for a protecting layer 38 (e.g. a glass layer) 
over the heating element, the index m for the thermographic material 5 and the index r for the platen 6 (e.g. a rubber 
roller or a drum). 

[0098] Now the step reading 'calculating transient temperature histories Th k ' in said thermographic material m for 
said sequences S k of time slices s, having different binary weights w, will be discussed in detail. 
[0099] The only time dependent quantity in a transient simulation of the thermal diffusion equation for the different 
slice combinations giving the temperature history in the thermographic material, is the specific heat production q(t) 
[expressed in J/m 3 ] in the nib material. 

[0100] We define a function u(t,t 8tart ,At) graphically depicted in Figure 8 as: 

Equation/: u(t, t^At) = 1 ift start < t < (g arf + At) 
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[0101] The heat generation in the heating element volume for a slice S k can then be expressed as: 



Equation 8 : 



V nib <=0 
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[expressed in W/m 3 ]. 

[0102] Herein P on [W] is the total power applied to a nib when a, equals 1 ; P off [W] the applied power to the heating 

element which is not actively steered (a,=0) but still might have some parasitic heat generation; and 

v n&> f m3 3 tne equivalent heating element material volume that will dissipate the power applied to it. 

[0103] The magnitudes of the applied powers P on and are known from the electrical configuration of the thermal 

head. 

[0104] In this way, the heat generation in the heating element is fully described by the slicing process. In total, one 
has 2 N available combinations for a certain fixed slice order. 

[0105] We have now an exactly described mathematical calculation scheme (FEM or FD) with known material prop- 
erties for every layer in our model (p, X., c), known boundary conditions and a known heat production in the nib layer 
q(t) for 2 N different combinations (see Fig. 14). The 1 dimensional heat diffusion equation (PDE) can now be solved 
for this problem: 



[0106] Details concerning the actual numerical details of these calculations are not given in this document, as this 
is a general mathematical issue. Many textbooks describe the theory about the FEM or the FD method that can be 
used for transient problems. Also, if necessary, other numerical techniques can be used. 

[0107] The solution will be available during the calculation in all nodes of the model. In practice, only the solution in 
the thermal sensitive layer of the thermographic material must be saved for later post processing. This should comprise 
at least one node, e.g. somewhere in the middle of the image forming emulsion. More nodes can be taken, depending 
on the evaluation function to be used later on. 

[0108] Now the step reading 'establishing a relationship between said transient temperature histories Th k and said 
output values D kcm will be disclosed thoroughly. 

In modelling the image forming behaviour in the thermographic material, special attention is focussed on the question 
what the pixel output is of the image forming process for a given transient temperature history. The theoretical function 
representing the pixel output will be denoted by Dh^ f(Th k ). Herein, at least two cases may be differentiated: 

» The nib is used to print pixel values, in which case the microscopic density (if necessary for a certain colour) in 

the pixel output is the representative target function for Dh k , or 
• The microscopic density is driven up to a saturating value, making only the size of the pixel formed by the nib being 

important (cf. raster imaging) and giving Dh k a dimensional value. Often only 1 dimension will be considered, e.g. 

the width of the pixel relative to the width of the nib. 

[0109] The input in the calculation of Dh k is the knowledge of the temperature history Th k in one or several points in 
the thermographic material. This relationship is in practice not obvious as it is based on chemical reactions and diffusion 
processes in the thermographic material. From the reaction kinetics and the optical properties of the chemical compo- 
nents, a temperature and time dependent relation with the pixel output Dh k can be established. 
[0110] We have found that a very good function for representing the pixel output Dh k is proportional to the maximum 
temperature reached in the thermographic material. In a thermographic system with the heat generated outside the 
thermographic material, the time constant of the heat present in some point, is generally greater than the time constant 
of the diffusion process governing the image forming. When a temperature T is reached, the temperature will be long 
enough around T for producing a pixel output which is independent in how that temperature T has been formed (within 
given constraints). 

[0111] In other cases, one might make a histogram of how long the thermographic material was exposed to a certain 
temperature. The longer it was exposed to temperature T, the more will be the contribution ofT. Some weighted integral 
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could then be used as the representing function of the pixel output Dh k . 

[0112] The above is illustrated with an example. For a given thermal head configuration, a transient simulation has 
been made for a slice with 8 weights having binary weights {1,2,4,8,16,32,64,128}. For the 255 levels, the maximum 
temperature Tmax occurring in the thermographic material has been calculated; see Figure 15, wherein 100 shows a 

5 maximum temperature relating to sequence weights W k , 102 shows a non-monotony in temperature at a sequence 
weight of 64, 104 shows a non-monotony in temperature at a sequence weight of 128, and ref. nr. 106 shows a non- 
monotony in temperature at a sequence weight of 1 96. This maximum temperature is in most cases representative for 
the pixel output Dh k on the thermographic material. One notices directly discontinuities around the major weight jumps, 
e.g. 63 to 64, 127 to 128, 195 to 196 and a lot of smaller jumps concerning the 15 to 16 lransition(the so-called "lowest 

10 nibble 01111 to 1000 transition' 1 ). 

[0113] The step of 'storing first entries for continuous and monotone output values Dh kcrT1 and second entries for 
corresponding sequences S k into a second table LUT2' now will be explained. 

[0114] As mentioned before, the pixel output D k is fully described by the printing sequence S k according to 

15 Equation 10: D k =F(S k ). 

From a software point of view, a sequence S k defines a certain slice order and a set of numbers defining active heating 
during a slice period. By selecting a set of sequences S k and ordering them into a table LUT, the index k of this table 
20 becomes a simple numerical unit that easily can be used in a software program. For the purpose of controllability, it is 
essential to make at least a continuous and monotone relationship between the index k of the lookup table and the 
corresponding pixel output obtained by activating the heating element by the sequence S k found at the table index k. 
[0115] This can be represented by a virtual function F cm , that is continuous as the function of the index k and mon- 
otone, according to: 

25 

Equation 11: if D k = F cm (k) and D f = F^I) and k<l then D k < D v 

[0116] The definition of F^ being continuous in the index k may to be formulated more specifically, as k is not a real 
30 number, but a natural number. Continuity of the output value D k is then formulated as giving for increasing k values, a 
continuous perception for the human eye, implying that for an increase in k of 1 , a maximum increase or decrease in 
the output D k is present, or: 

35 Equation 12: |D k -D k+1 | < AD max or ^{kyF^k+I^AD^, 

with AD max a maximum single step in output value. 

[0117] The function F^ has to be defined starting from a functional point of view and may comprise any shape, as 
far it fulfils Equation 11 and Equation 12 according to k; e.g. a linear shape, or parabolic, or hyperbolic, etc.... 
40 [0118] F cm and LUT2 are related to each other as: 

Equation 13: F^k) = F(LUT2[k]). 

45 Theoretically, the contents of the table LUT2 can be found by selecting a set of target values for D k and selecting for 
the second entry the sequence S k from which D k originates. But as mentioned before in the background of the invention , 
D k values are often unavailable as the energy used is not sufficient to produce a pixel output on the image forming 
material m. For the purpose of cross-talk or temperature offsets in the substrate or heats ink, LUT2 entries must be 
defined as well for these low energy cases. Whenever the energy is large enough to measure a pixel output, D k values 

so are mostly incorrect because of measuring errors, mutual cross-talk and unknown local substrate temperature. There- 
fore, we found that it is better to use the Dh k values. The step of storing first entries of Dh k values consists of selecting 
a continuous and monotone set GS' of Dh kcm values from the set GS of available and calculated Dh k values. From the 
relationship between Dh k and Th k , and Th k and S k , the second entry of the table can be found as being S k that will 
give a theoretical pixel output Dh kcm . 

55 [0119] In a further preferred embodiment according to the present invention, said step of storing into a second table 
LUT2 is replaced by a step of storing first entries for continuous, monotone and linear output values Dh kcm , and second 
entries for corresponding sequences S k in a third table LUT3, is explained by the following. 

[01 20] The mathematical process of compensation (cf . neighbour or cross-talk compensation , heatsink temperature 
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compensation) assumes that the selected sequence S k for heating element n produces a pixel output D k and further- 
more, that the principle of linearity is fulfilled or at least very good approximated. For the purpose of controllability, If a 
linear relation exists between the index k of the lookup table and the corresponding pixel output D k produced by the 
corresponding sequence S k of the second entry in that table, the principle of linear superposition can be applied when 

s describing the interaction of several sources contributing to the pixel output of a nib n. As an example, by cross-talk, 
the pixel n+1 when being printed, will give contribution to parasitic heat in pixel n. For a linear system, if pixel n+1 is 
printed using entry k of the table, a fraction f * k of the numerical value of k will be found in the pixel output of nib n. 
This can be compensated or when printing the pixel n by requesting an output found at the table index k - (f * k). The 
pixel output for pixel n will be the same as the output found for entry k when pixel n+1 is not printed. To realise this 

10 linear relationship, a table LUT3 will be constructed. The virtual function representing the image forming process F cm) 
which is continuous, monotone and linear according to table LUT3 conforms: 



D k = F cml {k)^aD k = F cml (a k\ae SRJ 

Equation 14: 

£>i = F cml {k x )and D 2 = /^,(* 2 )=> D x +D 2 = F cmJ {k { +* 2 ) 

20 with 

Equation 15: F cm iW = F(LUT3[k]). 

25 The function F cml must fulfil the continuity (Equation 11) and monotony constraint (Equation 12) and can be constructed 

as: 



Equation 16: F tem/ (*) = Dh^ *— + Dh of/set , 



with M the number of entries in LUT3, 
Dh max and Dh offsel a real constant. 
35 [0121] The first entry of the lookup table LUT3 is selected by selecting from the set GS a representative pixel output 
set GS" that contains Dh kcm , values according to the relationship: 

Equation 1 7: Dh kcmf = F cmf (k). 

40 

[0122] As a known relationship exists between the representative pixel output Dh kcm | and the temperature history 
Th k , and between Th k and S k , for a given first entry Dh kcm |, the corresponding sequence S k can readily be entered. 
[0123] There are two approaches available for constructing a linearised system of sequences. 
First of all, one could use a redundant system of slice weights. This means that the process of writing the time slice ts as: 

45 

Equation 18: t = J^a i w i t Scale 

/=o 

50 

is not unique, meaning that more than one set of {a,} exists in order to write t. One gets this situation e.g. when having 
two or more the same values of w, in the slice vector. The pixel output of the time t will be different for these cases, as 
the distribution over time of the weights will be different. A certain pixel output D k can then be realised in different 
55 manners. When selecting then a linear scale of Dh k , M in size, it will be no problem of finding a correct t value that can 
render this pixel size. M will be smaller than N. 

[0124] Another approach consists of keeping all the weights Wj different and searching a set of weights, in total M, 
being now equal to N, by simply changing the weight values until a monotonous and maybe even linear pixel output 
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Dh k is found with regard to the index k. For this, one can start form the binary weights -e.g. {1 ,2,4,8,1 6,32,64,...}- and 
for every transition from a set of smaller weights to one higher weight, trying to adopt this higher weight so that the 
transition happens smoothly and monotonous. Example: experiments show that it is better to start printing with the 
smallest weights and finish with the largest weights. Considerthe transition from {1,2,4,8,0, 0,0,...} to {0,0,0,0,16,0,0,...}. 
5 The slice set in the latter case will produce a pixel output with larger pixels than in the former case, with a perceptible 
jump, being not acceptable. A solution exists in dropping the value of the weight 1 6 until an acceptable increase of the 
size can be found. This process can be repeated for the higher slice weights. 

[0125] After having fulfilled this step, the lookup table LUT2 or LUT3 is constructed straightforward so that a contin- 
uous monotone output curve for a selected set of sequences S k is provided. 
10 [0126] In a further preferred embodiment of a method for thermal printing according to the present invention, the 
step of creating a time dependent power distribution P k = f(t) applicable to said heating elements H n , wherein said 
power distribution P k is defined by a known function of time t and having a total energy level determined by the parameter 
k, or Pj^fft.k), will be further established. 

[0127] In the state of the art thermal printing, using the heating element excitation layout of Fig. 4, the power applied 

is to the heating element has a time dependence as given by Equation 5. But as technology advances, more exotic control 
electronics will be available in the near future, possibly making the necessity of slicing obsolete. E.g. individual driving 
stages for each heating element could have a programmable timer, allowing to supply the heating element with P on 
during a programmable time. The heating element technology will also evolve. So, the use of an intense laser beam 
can be able to replace an heating element. I n that case, the time dependent power applied to the image forming material 

20 m, for a laser beam sweeping over it, may show a Gaussian time dependent profile and if necessary being repeated 
multiple times, giving a time dependent power applied to a particular pixel of the image forming material m, that is 
different from the constant voltage supplied heating element as in Equation 5. Therefore, the power supplied to a 
confined volume of the image forming material representing the pixel volume, can be denoted by a general time de- 
pendent function P k =f(t) [W], without having any constraints on the mathematical formulation of this function P k . As 

25 the problem of compensation will always present itself (increased ambient temperature or humidity, cross-talk between 
excited pixels etc.), the power delivered to the image forming material will be selectable and will mathematically be 
expressed by a parameter k, being a size limited integer number, having values 0 to K. For different values of k, P k 
will apply a different time dependent power to the image forming material and as a consequence, different values of 
the pixel output D k will result. The subject of the invention, as being explained and applied to a sequence S k of time 

30 slices S|, can readily be applied to P k . Using a calculation scheme, for the given material and geometrical set-up, the 
temperature histories Th k can be calculated in the image forming material m for all values of k. From these temperature 
histories Th k , a representative pixel output Dh k can be calculated, giving a set GS of Dh k values with their corresponding 
power excitation schemes P k . A subset GS' or GS" can be selected from GS as to give Dh k values that are continuous 
and monotone according to a predefined function F kcm , or continuous, monotone and linear, according to a linear 

35 relationship F kcml . 

(v) Illustration of the invention with a practical example. 

[0128] First, it is recapitulated that the core idea of the invention comprises: 

40 

• building a representative FEM or FD model of the thermal head, 

• making a transient simulation of the thermal diffusion equation for the different slice combinations, 

• running an evaluation function on the solution found, giving a representative output for the image forming process, 

• constructing an apt lookup table automatically so that one gets a continuous monotone output curve. 

45 

[0129] In the example of this section, the results depend of course very strongly on the type of thermal head used, 
and accordingly on the calculation scheme used for modelling the thermal head. But, in this context, it is just an illus- 
tration of the process of building an adequate lookup table. 

[0130] A Finite Difference FD approach has been used for a 1 D calculation scheme. A slice is used with 8 weights, 
so being {1 ,2,4,8,1 6,32,32,32}. The actual time slices are found by multiplication with a scale factort Scate , this for obtaining 
maximum density and pixelsize when printing a pixel with all slices on. The thermal head had such a configuration that 
when a pixel was driven during an on-slice, it dissipated a power P, and when not driven during a time slice, only P/9 
was dissipated (cf. Equation 5, P on and P off ). 

[0131] Several slices where taken with the same weight (here 32). This means that there are several possibilities to 
55 obtain a same summed weight, as an example: weight 64 can be made by {0,0,0,0,0,0,32,32} or {0,0,0,0,0,32,0,32} 
or {0,0,0,0,0,32,32,0}. But the pixel output will possibly differ, because the distribution of time the nib is dissipating heat 
is different. The choice of these sets of weights has been done deliberately to have "choices" when composing the 
table LUT2, representing the function F cm . Different weights will be available for making a same pixel, giving more 
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freedom in the construction of the lookup table. 

° In one experiment, a printout on the real printing device has been made (see Figure 16) with the binary weights 
from 125 {1 ,0,4,8,16,32,32,0} to 255 {1 ,2,4,8,16,32,32,32}. Herein, ref. nr. 110 is a pixel size relating to sequence 

5 weight, 1 1 2 is a discontinuity in pixel size at a sequence weight of 1 28, 1 1 4 is a discontinuity in pixel size and more 

particularly being an invisible pixel output zone, 116 is a continuous and monotone zone, 118 is a discontinuity in 
pixel size at a sequence weight of 196, and ref. nr. 119 is a saturation in pixel size at higher sequence weights. 
The pixel size is indirectly obtained from a macro density measurement. For the high excitation weights, a saturation 
occurs In the pixel output. This is due to a burning of the film, decreasing again the density of the film. There are 

10 also limits for a small nib to the absolute maximum pixel it can create. 

o In a second experiment, the numerical 1 D model was applied and as a function representing the pixel output, the 
peak temperature (T max ) occurring in the thermal sensitive layer, has been selected. The calculation result is graph- 
ically represented in Figure 17 showing calculated T max function for slice with weight vector {1 ,2,4,8,16,32,32,32}. 
15 Herein , following referral numbers are used: 120 is a partial curve of the maximum temperature in the thermo- 

graphic material m relating to sequence weights W K , 1 22 is a discontinuity in maximum temperature at a sequence 
weight of 128, 124 is a monotone zone showing some points of small discontinuities, 126 is a discontinuity in 
temperature at a sequence weight of 196, and 128 is another monotone and partly continuous zone. 

20 [0132] One finds a great similarity when comparing this picture (Figure 17) with the measured pixel size (Figure 16) 
if looking in the range from 125 to 255. Of course, the saturation effect 119 in the pixel output 110 for the high weight 
values is not visible in the T max calculation 120. 

[0133] Based on the values of T max , we can select a set of T max values which forms a linear monotonous vector of 
values. It does not make sense to try to form a vector with 255 elements. In our case, this is limited to 128 values in a 

25 first selection and after a few corrections limited to N=120. The selection is based on the principle to progressively 
select those sequences that give a continuous increase of the calculated T max value. Whenever a linear range is not 
available, one can take a range for some higher weight giving approximately the same T max value. Reference is made 
to Fig. 1 8 which shows the marked zones on the curve 1 32 relating T max to sequence weights and indicating the weights 
that have been selected for the LUT2 table. Herein, 130 represents a calculated Tmax for a weight system 

30 {1 ,2,4,8,16,32,32,32} with zones selected for the LUT2 table, 136 is a first selected zone of Tmax-curve , 137 is a 
second selected zone of Tmax-curve, 1 38 are third selected zones of Tmax-curve, and 1 39 are fourth selected zones 
of Tmax-curve. 

° In a third experiment as a final result, with the selected table LUT2, a printout has been made and the corresponding 
35 pixel size has been measured (Figure 20). Also the T max curve for the weights used has been shown in a diagram 

(Figure 19). Herein, ref. nr. 142 is a Tmax-curve related to LUT2-entries, and 144 is a pixelsize-curve related to 
LUT2-entries. 

[0134] One may remark that the pixel size does not behave fully linearly regarding the LUT2 index in Figure 20. Also, 
^0 the pixel sizes in this same figure are smaller than those in Figure 16. These effects are caused by an increased 
heatsink temperature, which makes a sole experimental attack of this problem practically infeasible as it is extremely 
difficult to run all the experiments under an isothermal condition. To compensate forthis change in heatsink temperature, 
one needs to install a heatsink compensation algorithm. 

45 (vi) Further preferred embodiments according to the present invention. 

[01 35] From the above, it may be clear that in a preferred embodiment of a method according to the present invention , 
said output values D K , and D kcm | relate to values of an optical density and/or to values of a pixel size to be 
reproduced on said thermographic material m. 
50 [0136] In a further preferred embodiment, said model describes a heat diffusion process in said thermal printing 
system. 

[0137] In a further preferred embodiment, said model is earned out by a heat diffusion partial differential equation 
PDE, comprising dimensions and thermal characteristics of all materials or layers contributing to said heat diffusion 
process. 

55 [0138] In a further preferred embodiment, said thermal characteristics comprise a specific mass p, a thermal con- 
ductivity X and a specific heat c. 

[0139] In a further preferred embodiment, said step of simulating said transient temperature histories Th k in said 
thermographic material m is carried out by solving said heat diffusion partial differential equation for ail sequences S k 



17 



EP 1 234 677 A1 



of time slices s { having different binary weights w ( . 

[0140] In a further preferred embodiment, said step of storing into said second table LUT2 comprises substeps of 

- making an explicit tabulation of different sequences S k and calculated output values Dh k , 

5 - selecting an ordered subset of output values Dh kcm showing a continuous and monotone relationship, 

storing into a second table LUT2 first entries for continuous and monotone output values Dh kcm and second entries 
for corresponding sequences S k . 

In still another preferred embodiment according to the present invention, a method for improving controllability 
of a thermal printing system comprising a thermal printer having a thermal head TH incorporating heating elements 

10 H n and a thermographic material m, said method comprising the steps of 

- creating a time dependent power distribution P k = f(t) applicable to said heating elements H n , wherein said power 
distribution P k is a known function of time t and having a total energy level determined by the parameter k, or P k =f 
(t,k); 

determining for every power distribution P k , a corresponding transient temperature history Th k = f(x,y,z,t) in said 
15 thermosensitive layer; 

. - calculating for each said transient temperature history Th k a corresponding output value Dh k , giving a set GS of 
output values Dh k with corresponding power distributions P k ; 

selecting from said set GS of output values Dh k , an ordered subset GS' of output values Dh kcni | realising a con- 
tinuous and monotone predefined output Dhk^F^fk); 
20 - storing into a table LUT2 first entries for said ordered subset GS' of output values Dh kcm and second entries for 
corresponding power distributions P k . 

[0141] In a further preferred embodiment, said step of storing into a second table LUT2 is replaced by a step of 
storing first entries for continuous, monotone and linear output values Dh kcm |=F cm ,(k) and second entries for corre- 
25 sponding power distributions P k in a third table LUT3. 

[0142] In a still further preferred embodiment, said output values Dh k , Dh kcm and Dh kcm , relate to values of an optical 
density and/or to values of a pixel size to be reproduced on said thermographic material m. 

(vii) Further applicability of a method according to the present invention 

30 

[0143] The method of the present invention is applicable for a wide variety of printing techniques. 

• In "Direct thermal printing", the method may be directed towards representing an image of a human body obtained 
during medical imaging and to a printing of medical image picture data received from a medical imaging device, 

35 e.g. a medical image camera. 

• Another application of the present invention comprises hardcopy printing for so-called non-destructive Testing 
(NDT), based on e.g. radiographic or on ultrasonic systems. Exemplary purposes of NDT comprise inspection or 
quality control of materials, welded joints or assemblies; development of manufacturing processes; experimenting 

40 in research; etc. 

• In another preferred embodiment of the present invention, the image data may be graphical image picture data 
received e.g. from a computerised publishing system. Further, a method according to the present invention also 
may be applied in graphic plotters, in chart recorders, in computer printers, etc. 

45 

• A method according to the present invention also may be applied in label printing by means of thermography. 
Thermographic label printing is described e.g. in WO 00/32403 (applicants Agfa-Gevaert and Esselte). More par- 
ticular applications within the field of label printing comprise cash register coupons, weight scale labels, etc. 

so • still another application of the present invention comprises thermographic printing of a bar code. Bar code tech- 
nology has become increasingly popular for identification and retrieval of all kinds of goods: from grocery store 
applications over medical laboratories, automobile industry, blood banks, electronics and libraries, etc. 

• Still another application of the present invention comprises thermographic printing of greeting cards and thermo- 
55 graphic printing of security documents. Security documents must be verifiable on their authenticity and comprise 

e.g. all kinds of identification documents such as passports, visas, identity cards, driver licenses, bankcards, credit 
cards, security entrance cards, and further value-documents such as banknotes, shares, bonds, certificates, 
cheques, gaming and lottery tickets and all kinds of entrance and transit tickets such as aeroplane tickets and 
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railroad tickets. 

» Thermal imaging according to the present invention can be used for production of both transparencies and reflec- 
tion-type prints. In the hard copy field, thermographic recording materials based on an opaque (e.g. white) base 
are used, whereas in the medical diagnostic field monochrome (e.g. black) images on a transparent base find wide 
application, since such prints can conveniently be viewed by means of a light box. 

«> A further application of the present invention provides means to assist blind people and persons with vision im- 
pairment. In "Research Disclosure" 32478 (dd. April 1991, page 268, disclosed anonymously) a method is de- 
scribed for obtaining embossed Braille-like output using special paper (e.g. coated with a layer of thermally ex- 
pandable plastic, e.g. a styrene compound) in a thermal printer. In a Braille coding, all characteristics (such as 
letters, numerals, punctuation marks, symbols, etc) are represented by specific combinations of six tangible points 
or dots. 

Heat supplied by the printhead of a thermal printer produces a perceptible output, because the localised heat 
produces a reaction in which the plastic layer creates raised dots. 

Alternatively, a computer program may be used to produce a series of raised dots representing data in graphical 
form, such as charts or maps. 

In order to get a continuous and monotone physical output for high resolution tactile graphics, a method according 
to the present invention may be used advantageously 

[0144] The present invention describes a technique of building a table that establishes a continuous and monotone 
relationship between a pixel output and an index value of a lookup table LUT. In a preferred embodiment, the table will 
be constructed as to give a more or less linear relationship, giving in fact an output system that has linear properties. 
This property of being linear is beneficial for any control algorithm and compensation technique that is used in controlling 
the pixel output. As an example, cross-talk between several heating elements (or nibs) presents a real problem in some 
commercial thick-film print heads. Cross-talk can be removed by performing a deconvolution technique on the image 
to be printed. Printing two pixels adjacent to each other in a same line will give then equal sized or equal dense pixels 
with regard to pixels being printed separately on different lines. 

[0145] For commercial thermal printers, it is often important to have a fast print rate. From a technological point of 
view, this requires that a nib can be heated in a short time, in order to produce a pixel output on the thermographic 
material. On the other hand, a rest heat has to be evacuated fast to a heatsink in order to cool down the nib at a 
temperature low enough to give no smog at a next line to be printed. Therefore, the heatsink may show an unwanted 
temperature rise, influencing directly a starting temperature of the heating element when printing a new pixel on the 
next line. It is important to have a model that correctly gives a relationship between the pixel output and an offset 
temperature being present in the heatsink, as to be able to compensate for this offset temperature. 
[0146] Further, it is important to indicate that for people skilled in the art, a so-called heating element may comprise 
e.g. a resistive heating element, an inductive heating element, a pyrotechnic heating element, or a high frequency 
heating element. 

[0147] Having described in detail preferred embodiments of the current invention, it will now be apparent to those 
skilled in the art that numerous modifications can be made therein without departing from the scope of the invention 
as defined in the appending claims. 



Parts list: 



[0148] 



1 

2 

3 

4 

5 

6 

7 

8 

10 

16 

18 

20 



a discontinuous curve 
a discontinuity 

a continuous, but non-monotone curve 
a non-monotony 

thermographic material m / thermographic imaging element le 
a platen 

hardcopy image 
slowscan direction 
thermal printer 
input data 
processing unit 
thermal print head TH 
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21 serial input line 

22 clock line 

23 latch enabling line 

24 strobe line 

5 25 shift register 

26 latch register 

27 AND-gate 

28 driver 

29 heating element 
10 30 control circuitry 

32 bonding layer 

34 substrate 

36 heating material 

38 protecting layer 

15 40 schematic for parallel to serial conversion 

42 schematic for serial outputs 

44 schematic for parallel inputs 

50 representation of function u 

52 a time slice having a weight w 0 and being in an "on-state" 

20 54 a time slice having a weight w, and being in an "on -state" 

56 a time slice having a weight w Nnibr ^ and being in an "on-state" 

58 a time slice having a weight w, and being in an "off -state" 

60 3D-representation of exemplary activation pulses 

61 activation of the heating elements during a first time slice 
25 62 second time slice activation of the heating elements 

63 last time slice activation of the heating elements 

64 first activation pulse for first heating element in a first slice 

65 first activation pulse for last heating element in a first slice 

66 last activation pulse for first heating element 
30 67 last activation pulse for last heating element 

68 activation pulse ('1') for nib N njbs -2 during third slice 

69 powerless activation pulse ('O') for nib 2 during third slice 
72 temperature for a single heat pulse of 2 mJ 

74 temperature for a first heat pulse of 1 mJ 

35 76 temperature for a second heat pulse of 2 mJ 

80 macrodensity relating to sequence weights 

82 jump in density at a sequence weight of 64 

84 jump in density at a sequence weight of 128 

86 jump in density at a sequence weight of 1 96 

40 88 saturation in density 

90 calculation scheme showing materials and properties 

100 maximum temperature relating to sequence weights {1 , 2, 4, 8, 16, 32, 64, 128} 

1 02 non-monotony in temperature at a sequence weight of 64 

1 04 non-monotony in temperature at a sequence weight of 1 28 

45 1 06 non-monotony in temperature at a sequence weight of 1 96 

1 1 0 pixel size relating to sequence weight 

1 1 2 discontinuity in pixel size at a sequence weight of 1 28 

1 1 4 invisible pixel output zone 

116 a continuous and monotone zone 

so 118 discontinuity in pixel size at a sequence weight of 1 96 

1 1 9 saturation in pixel size at higher sequence weights 

120 partial curve of temperature T max relating to sequence weights {1 ,2,4,8,1 6,32,32,32} 
122 discontinuity in maximum temperature at a sequence weight of 128 

124 a monotone zone with several points of discontinuity 

55 126 discontinuity in temperature T max at a sequence weight of 196 

1 28 another continuous and monotone zone 

130 calculated Tmax for a weight system {1 ,2,4,8,1 6,32,32,32} with zones selected for the LUT2 table 

132 T max -curve relating to sequence weights 
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1 36 first selected zone of T max -curve 

1 37 second selected zones of T max -curve 

1 38 third selected zones of T max -curve 

139 fourth selected zones of T max -curve 
5 142 T max -curve related to LUT2 -entries 

144 measured pixelsize-curve related to LUT2-entries 



Claims 

10 

1 . A method for printing an image using a thermal printing system comprising a thermal printer (1 0) having a thermal 
head TH (20) incorporating a plurality of energisable heating elements H n (29) and a thermographic material m 
(5), said method comprising the steps of 

is - supplying image data Id (1 6) corresponding to a plurality of output values D k to a processing unit (1 8) of said 

thermal printer, 

reading a first table LUT1 comprising first entries for desired output values D k and second entries for sequences 
S k of time slices s 1 of activation pulses, 

transformation of said image data Id into corresponding sequences S k of time slices of activation pulses, 
20 - providing said time slices of activation pulses to said heating elements H n of said thermal head TH, 

printing said image by transporting (8) said thermographic material past and adjacent to said thermal head 
and by activating (30) said heating elements H n of said thermal head, 

wherein said transformation of said image data Id comprises a step of establishing a relation between output values 
25 D kcm which are both continuous and monotone, and corresponding sequences S k of time slices Sj of activation 

pulses. 



2. The method according to claim 1 , wherein said step of establishing a relation between output values D kcm which 
are both continuous and monotone, and corresponding sequences S k of time slices Sj of activation pulses, com- 
30 prises the following substeps: 



building a calculation scheme (90) representative for said thermal printing system, said calculation scheme 
taking into account thermal characteristics of said thermal head TH (20), thermal characteristics of said ther- 
mographic material m (5), thermal characteristics of a platen (6) and specific heat productions q k as imposed 
35 by sequences S k of said time slices Sj having different binary weights w,; 

calculating transient temperature histories Th k in said thermographic material m for said sequences S k of time 
slices Sj having different binary weights w,; 

establishing a relationship between said transient temperature histories Th k and calculated output values Dh k ; 
storing into a second table LUT2 first entries for continuous and monotone output values Dh kcm and second 
40 entries for corresponding sequences S k . 

3. The method according to claim 2, wherein said step of storing into a second table LUT2 is replaced by a step of 
storing into a third table LUT3 first entries for continuous, monotone and linear output values Dh kcm) and second 
entries for corresponding sequences S k . 

45 

4. The method according to claim 1 , wherein said output values D k and relate to values of an optical density 
and/or to values of a pixel size to be reproduced on said thermographic material m. 



5. The method according to claim 2, wherein said calculation scheme (90) describes a heat diffusion process in said 
50 thermal printing system. 



6. The method according to claim 2, wherein said calculation scheme (90) is carried out by a heat diffusion partial 
differential equation, comprising dimensions and thermal characteristics of all materials or layers contributing to 
said heat diffusion process. 

55 

7. The method according to claim 2, wherein said thermal characteristics comprise a specific mass p, a thermal 
conductivity X and a specific heat c. 
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8. The method according to claim 6, wherein said step of calculating said transient temperature histories Th k in said 
thermographic material m is carried out by solving said heat diffusion partial differential equation for all sequences 
S k of time slices s, having different weights w h 

5 9. The method according to claim 2, wherein said step of storing into said second table LUT2 comprises substeps of 

making an explicit tabulation of different sequences S k and calculated output values Dh k , 
selecting an ordered subset of output values Dh kcm showing a continuous and monotone relationship according 
to a predefined function Dh kcm , 
10 - storing into a second table LUT2 first entries for continuous and monotone output values Dh kcm and second 

entries for corresponding sequences S k . 

10. A method for improving controllability of a thermal printing system comprising a thermal printer (10) having a 
thermal head TH (20) incorporating heating elements H n (29) and a thermographic material m (5), 

is said method comprising the steps of 

creating a time dependent power distribution P k applicable to said heating elements H n ; 
determining for every power distribution P k , a corresponding transient temperature history Th k in said thermo- 
graphic material m; 

20 - calculating for each said transient temperature history Th k a corresponding output value Dh k , giving a set GS 

of output values Dh k with corresponding power distributions P k ; 

selecting from said set GS of output values Dh k , an ordered subset GS' of output values Dh kcm realising a 
continuous and monotone output Dh^; 

storing into a table LUT2 first entries for said ordered subset GS' of output values Dh kcm and second entries 
25 for corresponding power distributions P k . 

11. The method according to claim 10, wherein said step of storing into a second table LUT2 is replaced by a step of 
storing into a third table LUT3 first entries for continuous, monotone and linear output values Dh kcm , and second 
entries for corresponding power distributions P k . 

30 
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